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Abstract 

The partition function for a canonical ensemble of 2D Coulomb charges in a 
background potential (the Dyson gas) is realized as a vacuum expectation value of 
a group-like element constructed in terms of free fermionic operators. This repre- 
sentation provides an explicit identification of the partition function with a tau- 
function of the 2D Toda lattice hierarchy. Its dispersionless (quasiclassical) limit 
yields the tau-function for analytic curves encoding the integrable structure of the 
inverse potential problem and parametric conformal maps. A similar fermionic real- 
ization of partition functions for grand canonical ensembles of 2D Coulomb charges 
in the presence of an ideal conductor is also suggested. Their representation as 
Fredholm determinants is given and their relation to integrable hierarchies, growth 
problems and conformal maps is discussed. 
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1 Introduction 



Statistical ensembles of 2D Coulomb particles (referred to also as logarithmic gases, (3- 
ensembles or Dyson gases) at a particular value of inverse temperature, = 2, are known 
to have remarkable integrable properties. Their partition functions, regarded as functions 
of properly chosen parameters of the background trapping potential (coupling constants), 
can be identified with tau-functions for hierarchies of nonlinear integrable equations such 
as 2D Toda lattice (2DTL) or Kadomtsev-Petviashvili (KP) hierarchies. Equivalently, 
the tau-function provides a generating series for correlation functions of the Dyson gas. 

The aim of this paper is to give an explicit representation of the partition functions 
as vacuum expectation values of certain operators constructed from free fermions, in the 
spirit of the Kyoto school [TJ [2] . Similar constructions in the context of random matrix 
models of different types were given in [31 HJ El E] . However, most of the previous studies 
were devoted to canonical ensembles with a fixed number of particles, N, which after 
F. Dyson [7j are customarily viewed as eigenvalues of a N x iV random matrix. The 
canonical partition function, Z^, is then tau-function of the 2DTL hierarchy with the 
discrete variable N. 

A natural question is whether this construction can be extended to grand canonical 
ensembles. Taking a weighted sum of the Zjv's, Z = J2n efiN %n, one obtains a quantity 
whose interpretation in terms of integrable hierarchies is presently not clear. A way to 
introduce a grand canonical ensemble of 2D Coulomb charges with transparent integrable 
properties was first suggested in [8J. The idea is to consider the logarithmic gas with 
varying number of particles in the presence of an ideal conductor, so that each particle 
interacts not only with other particles of the gas but also with their "mirror images" 
of opposite charge. The grand canonical partition functions of such systems appear to 
be tau-functions of the 2DTL or KP hierarchies (depending on whether the conductor 
fills a disk or a half-plane), whose "times" again serve as parameters of the background 
potential but in a different manner. The fermionic operator construction of these tau- 
functions is technically even simpler than that for canonical ensembles. 

Our motivation comes from growth problems of Laplacian type such as viscous flows in 
the Hele-Show cell [9], [10]. At zero surface tension, the model has an integrable structure 
of the 2DTL hierarchy in the zero dispersion limit [llj. From mathematical point of view, 
the same integrable hierarchy stays behind some classical problems of complex analysis: 
parametric deformations of conformal maps [T2], boundary value problems for Laplace 
operator [131 E] an d the inverse potential problem in two dimensions [15J. Switching 
on the dispersion thus coming back to the original 2DTL hierarchy, one obtains a non- 
trivial deformation of this integrable structure. However, the corresponding deformation 
of Laplacian growth and, what would be even more intriguing, of the above mentioned 
problems of complex analysis, is not yet formulated explicitly in proper terms. 

The theory of logarithmic gases provides another view on these matters [TBI E] • In 
the thermodynamic limit, when the number of particles is very large, the logarithmic gas 
in the leading approximation macroscopically looks like a charged fluid with continuous 
density. The equilibrium state of the system is a result of competition between the 
mutual repelling of particles, which tends to remove them to infinity, and the external 
force which attracts each particle to local minima of the background potential. Typically, 
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in the equilibrium the fluid occupies some compact domains in the plane around local 
minima of the background potential, which we call droplets. The shape of the droplets 
is determined by their total charge and by the profile of the trapping potential. When 
one increases the total charge keeping the potential fixed, the droplet changes its shape 
according to the growth law specific for the Laplacian growth processes (the Darcy's law). 
At the same time, the integrable deformation of the N = oo picture is naturally build 
in the 2D Coulomb gas. We believe that the reformulation of the model in terms of free 
fermions may help to clarify the geometrical meaning of this deformation. 

Another interesting question is what kind of growth problems emerge in the thermo- 
dynamic limit of the grand canonical ensembles. We give a partial answer restricting 
ourselves to the grand canonical ensemble in the upper half-plane, in which case the 
growth problem is equivalent to Laplacian growth of "fat slits" considered in [T8] . 

2 Free fermions and tau-functions 

In this section we recall the free fermionic construction [TJ [2] of KP and 2DTL tau- 
functions. 

Let ip n , t/)*, n e Z, be free fermionic operators with usual anticommutation relations 
[ip n ,ipm]+ = [i>nii ) m]+ = °> ['0n,V'm]+ = <W They generate an infinite dimensional 
Clifford algebra. We also use their Fourier transforms 

ij>(z) = £ ^ k z\ r{z) = E ^ (2.1) 
fcez fcez 

which are regarded as free fermionic fields in the complex plane of the variable z. 

Neutral bilinear combinations ^mn'4 ) m'4 , n °f the fermions generate the Lie algebra 
g/(oo). Normal ordering of the generators (see below) allows one to consider certain infi- 
nite sums as well. Exponentiating these expressions, one obtains an infinite dimensional 
group which is a central extension of GL^. Let us call elements of the Clifford algebra of 
the form g = exp (J2mn A mn ip m ipn) group-like elements. A characteristic property of the 
group-like elements g is that g^ n 9~ x is a linear combination of ipj's and similarly for ?/>*: 
QTpnQ' 1 = Hi^iRim gi'nd 1 = ^i^i^nii wnere the matrix R is determined by the ma- 
trix A. Of particular importance are the group-like elements obtained by exponentiating 
the operators 

J+ — ^ tkJkt J- — t-kJ-k (2-2) 
k>\ k>l 

where 

J k = E Wj+k 

are Fourier modes of the "current operator" and t k are parameters (called times). In 
general, they are complex numbers. It is convenient to denote the collection of times 
with positive (negative) indices by t + = {ti, t 2 , ■ • ■} and £_ = , t_ 2 , • • •} respectively 
and to introduce their generating function 

iit ±1 z) = Y,t±kZ k . 

k>\ 
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Figure 1: The Dirac sea |0). Filled states are shown in black, empty ones in white. 

It is easy to check that the fields ip(z), ip*{z) transform diagonally under the adjoint 
action of the group-like elements e J+ , e J ~ : 



e J± ijj(z)e- J± = e ?(t± - 2±1 V(z) 
e J± ip*{z)e- J± = e-^' z±1 ^*{z). 



(2.3) 



In terms of polynomials Pk{t±) defined by 



e «'±.*) = J2Pk(t±)z k (2.4) 

fc>0 



the corresponding formulas for if) n ,ifj^ can be written as 



k>0 



(2.5) 



e j ±r n e- J± = Y,r n±k Pk(-t±) 

k>0 



Next, we introduce a vacuum state |0) which is a "Dirac sea" where all negative mode 
states are empty and all positive ones are occupied (Fig. CQ: 

V>«|0) = 0, n<0; <|0) = 0, n > 0. 

(For brevity, we call indices n > positive.) With respect to this vacuum, the operators 
ip n with n < and ■?/>* with n > are annihilation operators while the operators ^* with 
n < and ip n with n > are creation operators. Similarly, the dual vacuum state has 
the properties 

(0|< = 0, n<0; (0|^n = 0, n > 0. 

We also need "shifted" Dirac vacua \n) and (n| defined as 



\n) 



(n 



V>»-i • • ■V'lV'o |0) , ^>0 

C---C 2 Ci|0), n<0 

(0|#i ^>0 

(O|'0-l'0-2---V'n ) ™<0 



The vacuum expectation value (0| . . . |0) is a hermitian linear form on the Clif- 
ford algebra defined on bilinear combinations of fermions by the properties (0 1 0) = 1, 
(0| ipntpm |0) = (0| |0) = for all m, n and 

(0\M^\0) = 5 mn form<0, (0|^ n ^|0) = for m > 0. 
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Its extension to the whole algebra is given by the Wick's theorem: Let Wi be arbitrary 
linear combinations of ipj and ip*, then (0| W\ . . . w 2n +i 1 0) = and 

(0| Wl . . . W 2n |0) = £(-l) P(<T) (0| W a{1) W a{2) |0> <0j W a( 2n-l)W*(2n) |0) . 

a 

Here the sum is over permutations a of 1, . . . , 2n such that a(l) < <r(2), <r(3) < cr(4), 
. . ., a(2n — 1) < <j(2n) and a{\) < cr(3) < . . . < a(2n — 1) and P(a) is the parity of a. 
The Wick's theorem is often used in the following form. Let Wi be linear combinations 
of ip/s only and w* be linear combinations of ip*'s only, then 

(0\w 1 ...w n w* n ...w* 1 \0) = det (0\wiW*\0). 

i,j=l,...,n 

Let us also give an explicit formula for the expectation value of products of the fields 



n t 

(N\ i>( Zl ) . . . i>(z n )r((n) ■ ■ ■ \N) = U(zi/(i) N ■ det • s 



i,j Zi - (j 



TT ^iSl TT ( Zi Z j)(Cj Ci) 

Mi:; Q) \} 3 ( Zi - Q( Zj - QY 



(2.6) 



One may define the normal ordering *(...)* with respect to the Dirac vacuum |0) when 
all annihilation operators are moved to the right and all creation operators are moved to 
the left taking into account their anticommutativity under mutual permutations. For ex- 
ample, = -ipii^l = ipi^i-l, and, more generally, \^* m ^ n \ = ^* m ^ n -{^\^* m ^n |0). 
We also note the identities 

e aM* k = 1 _|_ ( e a _ l)^* k = ' e (e a -l)i> k r k • for k > 0, 

(2.7) 

e <*l>Wk = 1 _|_ ( e a _ l)^*^ fe = . e (e«-l)^^ • for k < 0. 

As it has been established in the works of the Kyoto school, the expectation values 
of group-like elements are tau-functions of integrable hierarchies of nonlinear differential 
equations. In particular, 

r n (t+,f_) = (n\e J+ ge' J - \n) (2.8) 

is the tau-function of the 2DTL hierarchy meaning that it obeys the infinite set of Hirota 
bilinear equations for the 2DTL hierarchy the simplest of which is 

dtiTndt^Tn - T n d tl d t _{T n = T n +\T n -\- 

In a similar manner, 

r(t + ) = (0\e J +g\0) (2.9) 

is the tau-function of the KP hierarchy. 

It is worthwhile to note that the group-like elements of the Clifford algebra can be 
also written as normal ordered exponents of the form 

s = :exp(£i m „VmC) : (2-io) 



For regular (invertible) elements g, this normal ordered expression is just equal, modulo 
a constant factor, to an exponent of a fermionic bilinear form without normal ordering 
but with some other matrix A (an example is provided by (12.71) . see also fTS] for a more 
general case). If g is invertible, then the tau-function at any n is not identically zero. 
However, there is an important class of tau-functions which vanish identically at some 
values of n. They correspond to singular elements g which are not invertible (and thus 
do not belong to a group in a strict sense) but still can be represented in the normal 
ordered form f 1 2 . 1 1) . 

Here are two important examples of such singular elements: 

p + = :ex P fe^*) : = IK 1 - = UMl 

\i<0 } i<0 i<0 

(2.11) 

p_ = : ex P f -2> 4 < ) : = na - m:) = n ^- 

\ i>0 / i>0 i>0 

In a sense, these operators are projectors to positive and negative modes respectively. 
Their properties (extensively used in what follows) can be easily seen from the definition. 
The both operators obey the projector property: P± = P±. The operator P + kills 
negative creation modes standing to the right and negative annihilation modes standing 
to the left and commutes with all positive modes: 

p+r k = ^ k p+ = o, k < o, 

(2.12) 

[P + ,^] = [P+,^.]=0, k>0. 

The operator P_ kills positive creation modes standing to the right and positive annihi- 
lation modes standing to the left and commutes with all negative modes: 

P^ k = r k P- = 0, k > 0, 

(2.13) 

[P_,^] = [P+,Vfc] = 0, k<0. 

From this it is obvious that P+\n) = at n < and P+\n) = \n) at n > 0. Similarly, 
P-\n) = at n > and P-|n) = \n) at n < 0. Somewhat less obvious properties (also 
used in what follows) are P + e~ J ~ |0) = |0), (0| e J +P + = (0|. 

As an example, let us calculate the tau-function corresponding to the singular element 
P + : Tn = (N\e J +P + e- J - \N). Here we follow [3]. First of all, it is not difficult to see 
that tn = at N < and tq — 1. For iV > 1 we have: 

rjv = (0| r N -i ■■■^o e J+ P + e- J -ip • • • Vtf-i |0) • 

To proceed, it is convenient to use the short hand notation ip n {H) = e H ip n e~ H , ipn(H) = 
e H ipn e ~ H > where H is any operator from the Clifford algebra, then we can write 

r N = (0| e J +r N -i(-J+) ■ ..r (-J+)P+P+M-J-) ■ ..ipN-i(-J-)e- J - |0> . 
Equations (12.51) imply 

k>0 

(2.14) 

^n(- J +) = Y.^*n+kPk{t+). 
k>0 
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so the operators ip n (—J^), ip^(~J + ) in the formula for t n contain only positive modes 
and, therefore, commute with P + . Moving one P + to the right and another one to the 
left, and using the properties mentioned above, we obtain 

r N = {o\r N -i(-j + )---ro(-j + )M-J-)--^N-i(-j-)\o) 

= det (0|^_ 1 (-J + )V fe -i(-^-)|0) 

by the Wick's theorem. The expectation value under the determinant can be represented 
as a contour integral as follows: 

(o\rj(-j+)M-J-) io> = E va(t + ) Pb {-t-) (oiv>; +( a +& |o> 

a,b>0 

= E Pa(t + ) Ph (-t-)5 J+a , k+h = I Z ^ e m,)-^)J±_ 

a,b>0 J \ z \ =1 2nlZ 

The whole determinant can then be written as an iV-fold contour integral: 

TN = ^j...f I] " - ^) I! ^'^-'^^ ■ ( 2 - 15 ) 

When t^k = —tk, the expression £(t+,z) — £(£_, l/z) is purely real for z on the unit circle 
and tn coincides with the partition function of the unitary random matrix model written 
in terms of the eigenvalues. In this form, it can be treated also as the partition function 
of a canonical ensemble of iV 2D Coulomb particles confined on a circle. 



3 Partition function of the 2D Coulomb gas as a tau- 
function: canonical ensemble 

Let us fix an arbitrary measure d/i(z) in the complex plane and consider the following 
group- like element: 

g = : exp ( [ J, + (zW + (l/z)dp(z) ; (3.1) 

Here 4>+(z) = E i J nZ n , i>* + {z) = E "^n- 2 "™ are truncated Fourier series containing only 

n>0 n>0 

positive modes. Obviously, g commutes with P + . Expending the exponent into a series, 
one can represent go in a more explicit form: 

00 1 r 

go = E — 1 / m ^+( Z ^- ■ - , ip+{z m )P-'ip* + (l/z m ) . . .^+(l/^i)d/Xi. ..dfJLm, (3.2) 



where dfij = d/i(zj). 

Let us consider the expectation value 

r N (t + , t_) = (N\ e J +g P + e- J - \N) (3.3) 
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and apply to it a chain of transformations similar to the ones used in the simpler case 
go = 1. Again, tjv = at N < and r = 1. For iV > 1 we have: 

r N = (0|^_ 1 ...^e J +^ P + e- J -^ ...^_i|0) 

= (0| e J +r N -i(-J+) ■ ..^{-J + )P+goP+M-J-) ■ .^N-i(-J-)e- J - |0> 

= <0| ^-i(-^t-) • • • r o (-J+)9oM-J-) ■ ■ ■ ^n-i(-J-) |0> . 
Substituting the explicit form of go, we get: 

t n = —1 I m d ^- ■ ■dfi m (0\ip* N -i(-J+) ■ ■ .^o(- J +)^+(^i) • • -Mzm) 

m>0 

X P^* + {l/z m ) . . . ■ ■ ■ TpN-l(-J-) |0>. 

The next step is to notice that only the term with m = N contributes to the sum and 
all other terms vanish. Indeed, at in > N the state 

V>;(l/^m) • • ■ rl>*+0-/zi)M-J-) ■ ■ ■ 4>n-i(-J-) |0> 

is in fact the null state because the number of annihilation operators exceeds the number 
of creation operators while at m < N the operator 

p_^;(i/z m ) . . . r+(i/h)M-J-) ■ ■ ■ iI>n-i{-J-) 

is in fact the null operator because P_ multiplied by the uncompensated positive ?/>-modes 
from the right gives (see ( 12 . 1 3[) ) . Therefore, the expression simplifies to 

™ = ivT Ic N dfil " ' dfiN ^* N -^~ J+ ^ " " " ^o(- J +)^+(^i) • • • M*n) 

x p_^;(i/%) . . . r + (i/h)M-J-) ■ ■ ■ ^n-i(-j-) |o). 

Since there are as many annihilation operators to the right of P_ as creation ones, the 
state that they produce from the vacuum is proportional to the vacuum state itself, i.e., 

v;(i/%) . ..r + (i/h)M-J-) ■ ..^n-i(-j-) io> = |o> c N , 

where the constant Cn is 

c N = (o\r+('i-/z N )...r + ( i /^)M-j-)--^N~i(-j-)\o) 

= det (0|^;(l/%)V fc -i(-^-)|0). 

(0|V;(l/^fc-l(-J-)|0)= £ ^Pa(-t-)(0\^k+a-l\0) 

a,l>0 



a>0 
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the constant Cn is explicitly given by 



N 

C N = A N ( Zi )l[e-^->*'\ (3.4) 
1=1 



where we use the convenient short-hand notation for the Vandermonde determinant: 

M^) = 1 det 2 ,(4- 1 )=II^- 



Now, it remains to calculate 

<0| Pn-i(-J+) ■ ■ ■ r (-J+)Mzi) ■ ■ ■ M*n) |0> = det (0| ^(- |0> 

N 

= det zl l e^ Zk) = A w (^)TTe €( ' + ' 2i) , 

Kj,k<N K , , 

which can be done in a completely similar manner. Collecting everything together, we 
obtain the result: 



1 r N 
r N (t + ,t.) = — IAnWI'U^"™*-'*^). (3.5) 

iV . J<L l=1 



Assume that djj,{z) = e~ u( - z ' z ^d 2 z is a smooth measure on the plane and t_k = —i^, 
then the expression £(£+, z) —£(£_, f) is purely real and the integral (13. 5p has a physical 
interpretation as the partition function of a canonical ensemble of iV identical Coulomb 
particles in the plane in the background potential W(z,z) = —U(z,z) + TR,e^ k t k z k : 



1 r N 



It is proportional to the partition function of the ensemble of normal random N x N 
matrices $ [20]: Zn oc J D§ e trw(*,*t)^ w ^ z . being their eigenvalues. 

If the measure d\x is concentrated on a curve T C C, then the 2D integrals J c (. . .)g? 2 z 
are reduced to ID integrals / r (. . .)\dz\ along T. This means that the 2D Coulomb particles 
are confined to the curve T. For particular choices of T the integral (I3.5P yields the 
partition functions of random matrix models of certain types in terms of eigenvalues. 
For example, if F is the real line, one obtains the partition function of hermitian random 
matrices and if T is the unit circle, the integral (I3.5P becomes identical to (I2.15P which 
is the partition function of unitary random matrices. 

Let us show that for axially symmetric measures dfi the operator representation (I3.3P 
is equivalent to the one suggested by A. Orlov et al [21]. For an axially symmetric 
measure, the bilinear form in the fermion operators in (13. ip becomes diagonal: 

$+{z)il>X(l/z)dii{z)= £ 2Ye-^¥2 = £M^ 

m,n>0 C n>0 

where 

h n = I \z\ 2n e- u ^d 2 z, 



9 



(we assume that the measure is smooth with U(z, z) = U(\z\)), so in this case 

9o = I exp 5Z(/i n -l)Vv0* K = exp ^ \ogh n ip n ip* (3.7) 

\n>0 / \n>0 / 

and the tau-function (I3.3P does have the form (N\ e J+ e x e~ J ~ \N) with X =y^ Xj * ipjip* * 

studied in [21]. More precisely, it corresponds to a singular limit of the latter with 
Xj = loghj for j > and X, — >■ +oo for all j < 0. Indeed, writing 

e x = U(l + (e x * - 1)^*) ■ (l + (e"*' - 
i>o i<o 

we see that the first product is equal to n (1 + (hj — l)ipjip*) = g while the limit of the 

j>0 

second one is the singular operator P + . 

An important example is U(z, z) = c\z\ 2 , then h n = f[c~ n ~ l n\ = Trc~ n ~ 1 r(n + 1) and 
X n = — (n + 1) logc + logT(n + 1) (the common constant log7r is irrelevant because the 
operator ^ '.ipjipjl commutes with all elements of the Clifford algebra). Note that the 

analytic continuation of this formula to negative values of n with the help of the gamma- 
function automatically implies the required singular limit X n = +oo at n < 0. The 
tau-function (13. 5ft for this case has the following expansion in Schur functions (21} [22] : 

N 

r N (t + ,t.) = rc N c- N ^/ 2 1] r(*) ■ J2c-^(N) x s x (t + )s x (-t-). (3.8) 

fc=l A 

Here A denotes the Young diagram with £(A) rows of lengths A x > A 2 > . . . > A^m > 0, 
| A | is the total number of boxes in A, 

/(A) 

(N) x := (N + 1 - i)(N + 2 - i) . . . (N + X t - i) 
i=i 

and the sum runs over all Young diagrams including the empty one (its contribution 
is 1). The Schur functions s x are defined in the standard way as determinants of the 
polynomials Pk(t) ( 12. 4ft : 

s x(t) = deb p Xi -i +j {t). 

1<«,J<£(A) 

The dispersionless limit of the tau-function (13. 5p is achieved via rescaling of times 
tk — Tk/h, N = To/h and tending h — > 0. Clearly, this implies N — > oo. However, in 
order to get a meaningful limit the measure d\x should be chosen appropriately. In the 
case of a smooth measure one should set dfi(z) = e~K U( - z ' z ^d 2 z, then 

r N (t + , £_) = exp ^-J + Ofi- 1 )^ , h -> 0. (3.9) 

The function F — F (. . . ,T_i,T ,Ti, . . .) is what is called "dispersionless tau-function" 
[23| |24"] . When the reality condition T_fc = — Tk is imposed, this function encodes a formal 
solution to the inverse potential problem in 2D and admits a nice geometric/electrostatic 
description [T5]. Set a(z, z) = ^d z d s U(z, z) to be density of the background charges. We 

10 



assume that a > 0. Given a real positive T and complex TVs, k > 1, it is a subject 
of the inverse potential problem to find a domain D in the complex plane such that T 
is the total charge contained in D and TVs are harmonic moments of its exterior with 
respect to the density a: 

T = [ cr(z,z)d 2 z, T k = -\( z- k a(z,z)d 2 z, k > 1. (3.10) 

JD k JC\D 

Leaving aside the very difficult questions about existence and uniqueness of the solution, 
we simply assume, just for an illustrative purpose, that we are in a situation when D is 
a compact connected domain (containing the origin). Then F is given by 



F = - a(z, z) log z- 1 - C 1 *((, d 2 zd 2 C (3.11) 

JDJD 

which is basically the electrostatic energy of D filled by electric charge with density a 
(and with a point-like charge at the origin). In the case U(z,z) = zz, a = 1/tt one 
obtains the "tau-function of analytic curves" [25] which encodes the integrable structure 
of parametric families of conformal maps [12] and the Dirichlet boundary value problem 
in C \ D [131 CH]- The analytic continuation of the function Fq to general values of 
(not necessarily constrained by the condition = — T*.) also has a geometric meaning 
in terms of pairs of conformal maps [26J. As the results of [27] suggest, higher terms 
in the ^-expansion of the tau-function (13.91) may be related to spectral invariants of the 
domain D. 

For tau-functions of the form (13.51) with a singular measure dfi concentrated on a 
contour, the dispersionless limit does not require introducing any ^-dependence of the 
measure. The leading h — > behavior has the same form (13. 9ft and there is an integral 
representation for F similar to (13. lip with the domain D being replaced by a segment 
of the curve. However, the parameters 7* in this case do not admit a direct geometric 
interpretation. 



4 Partition function of the 2D Coulomb gas as a tau- 
function: grand canonical ensemble 

A direct attempt of passing to a grand canonical ensemble via Z — ^ e^ N with 

N>0 

given by (I3.6P leads to a quantity whose interpretation as a tau-function of an integrable 
hierarchy is not known. A way to introduce a grand canonical ensemble with good inte- 
grable properties, first suggested by I. Loutsenko et al [8], is to consider the logarithmic 
gas in the presence of an ideal conductor, so that each particle interacts not only with 
other particles of the gas but also with their "mirror images" of opposite charge as well 
as with its one image. For the reflection to be globally well-defined the conductor should 
be a disk. 

Let D be the unit disk and O* its exterior. We assume that particles of the gas with 
complex coordinates z\ occupy D* and © is an ideal conductor. Then the mirror images 
are inside D at the points l/z{ and the Coulomb energy of the system is 

N N 
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The grand canonical partition function of the system of identical particles in a background 
trapping potential W (which is supposed to include the chemical potential) is defined as 



Z =f:^ [ ■■■[ e E " +w "d 2 Zl ...d 2 z N 

N=0 iv ! Jo* Jo* 

In a more explicit form, it reads 
Z 



N 



oo i 

E 



7V=0 



N 

n 

j<k 



Zj Zfc 



ZjZk — 1 



2 N 

n 



W N = Y,W{z u zi) 
1=1 



Zl \ e W(zi,zi) 

d z h 



N 2 -i 



(4.1) 



In order to ensure convergency of the integrals, the function e w( - z ' z ^ should vanish at 
|*| = 1 as (\z\ 2 - l) a with a > 0. 

Let us fix an arbitrary measure dfi(z) in D* and consider the following group-like 
element: 



G 



cxp 



(4.2) 



Here we use another normal ordering * (. . .) * , the one with respect to the completely 
filled vacuum | — oo), which means that the ?/>-modes are moved to the left while ^*-modes 
are moved to the right. Expanding the exponent into a series, one can represent G in a 
more explicit form: 



N=0 



1p{zi) . . .1p(z N )lp*(l/z N ) . ..^(l/z^dfMi . ..d/Ir, 



(4.3) 



where dfij = dfi,(zj). Using (I2.6P it is straightforward to see that the expectation value 



r^ G) (t + ,t_) = (n\e J +Ge- J - \n) 
has the structure of partition function of a grand canonical ensemble: 



where 

u(t] z, z) 



■ J2k>i ktkt - 



E 

N>0 



N 

n 

j<k 



Zj Zfc 



Zj z k — 1 



2 N 

n« 

i=i 



u(t;zi,zi) 



dfi(zi) 

N 2 -i' 



(4.4) 



(4.5) 



to log \z\ 2 + £(t+, z) + £(t_, l/z) - £(t+, l/z) - £(t_, z) 



t \og\z\ 2 + Y.( t k( 



z k - z~ k ) 



tJz k - z- k \ 



k>l 



(4.6) 



and n = t . Clearly, at = — u)(t\ z, z) is real and can be interpreted as a harmonic 
part of the background potential in D*. More precisely, the identification with (14.11) goes 
as follows: set W(z,z) = —U(z,z) +u(t;z,z) — log|z|, then dfx(z) = e~ u ( z ' z ' ) d 2 z. It is 
convenient to redefine the tau-function by extracting the simple factor in front of the 
sum in ?P(t + ,t_) = e^ 1 HfeLfi rf^ + , then Z = (*+,*_). 

Note that the tau-function is formally a oo-soliton tau-function with momenta of 
solitons distributed with the measure d\x. It is the Fredholm determinant of an integral 
operator: 

fP(t + ,t_) = det(l + K). (4.7) 
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Here 1 is the identity operator and the operator K acts to functions on D* as follows: 

Kf(z)= / £W—-dtiC). 
Jo* zQ — 1 

To see this, we again calculate the expectation value of ( 14. 3 I) with the help of (12.61) but 
now using the determinant representation of each term. This gives an expansion of the 
Fredholm determinant. 

As an example, consider the case when the measure dfj, is concentrated on a circle of 
radius e e , e > 0. Set z = e e+l< ^, then d^(z) = e e d(p and the sum in the r.h.s. of (14. 5 j) 
becomes the grand canonical partition function of the charged particles on the circle 
\z\ = e e in the presence of the ideal conductor: 



sm 



sin fr-y- 2 * 6 



2 TV 



i=i 



where the chemical potential fi is expressed through to as /i = 2eto — log (e e — e~ e ) and 

V(0) = 2 sinh ( efc ) ^ e ^- 

This partition function is a compactified version of the one considered by V. Kazakov et 
al in [28]. As a function of the "times" it is tau-function of the 2DTL hierarchy. 

A slightly more general setting is to take the disk-like conductor of an arbitrary radius 
R and to consider the grand canonical ensemble of charged particles in its exterior. In 
fact for any finite R this hardly brings anything new because the systems at different 
R can be transformed into each other by a simple rescaling of z and by redefining the 
background potential. However, in the singular limit R — > oo a new grand canonical 
ensemble emerges, which can be introduced independently without a reference to any 
limiting process. It is defined in the upper half plane H, with the lower half plane 
being an ideal conductor. Its partition function is tau-function of the KP hierarchy. 
The limiting procedure is rather sophisticated and will not be described here. We only 
remark that it is similar to the one developed by E. Antonov et al [29] for the transition 
2DTL — > KP in the continuum limit of the 2DTL hierarchy. 

Let us give the main formulas related to the grand canonical ensemble in EI and its 
operator realization. Fix an arbitrary measure d^{z) in HI and consider the following 
group-like element: 

G n = I exp (ijjj(z)iP*(z)z- l dfM m (z)j *. (4.9) 

The expectation value 

r(t) = (0|e J+ G H |0) (4.10) 

is KP tau-function depending on the "times" t = {tx, . . .}. (For the abuse of notation 
we use the same letters as in (14.41) but one should remember that their meaning is 
different, see below.) A calculation similar to the one for D* yields 



-w=E^/ H -y H n 

Af>0 iV • Ja JML j<k 



N _ „ 2 N 



Zj Zfc 



Zj Zfc 



TT e ««;«<)-e(«;3) d ^( z i) u n \ 

ZLmzi ' ' 



i=i 
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a) 



b) 



Figure 2: A fat slit B a) in the upper half plane, b) in the exterior of the unit disk. 



which looks like a oo-soliton tau-function with continuously distributed momenta. (The 
analogy with soliton solutions of the KP hierarchy was first noticed in [8].) In order to 
ensure convergency on the real line, the measure should vanish there as dfi^_(z) cx (Imz) a 
with a > 0. Assuming that all the t^s are purely imaginary and dfxu(z) = e~ u ( z ' z >cPz, 
this tau-function is identical to the partition function of the grand canonical ensemble of 
2D Coulomb particles in EI with the background potential W(z, z) = —U(z, z) +£(£, z) — 
£(t,z). Choosing the measure dfim concentrated on the line x + is, one obtains from 
(14. lip a non-compact analog of (14.81) . which is precisely the partition function considered 



m 



28 . 



The dispersionless limit of the tau-function (14. lip is again achieved via rescaling of 
times t k = T k /h, setting the measure to be dfj,u{z) = e~^ u ^ z ' z ^d 2 z and tending h — > 0. 
In contrast to the canonical case, the mean total charge of the system can not be taken 
arbitrary but is fixed by the equilibrium condition. In other words, at h — > the leading 
contribution to the sum ( 14. lift comes from its maximal term. One has: 



r(t) = exp 




h -»■ o. 



(4.12) 



The function Fq = Fq(Ti,T2, . . .) is the tau-function of the dispersionless KP hierarchy 
[23~t 124"] . When the reality condition TZeT k = is imposed, this function admits a nice 
geometric/electrostatic description given in [18] for the particular case o = 1/tt. Let 
(t(z,z) = ^d z dzll(z, z) be density of the background charges, as before. Given purely 
imaginary Tjt's, k > 1, a version of the inverse potential problem in the upper half plane 
allows one to find a domain B (a "fat slit" [T8], see Fig. |2£i) in H such that 



T\ = —2iXm / a(z, z)d 2 z, 
Jb 



Tk = —Ira 



H\B 



z- k a(z,z)d 2 z, k>2. 



(4.13) 



Then F is given by 



BJB 



a(z, z) log 



C 



^(CO^C, 



(4.14) 



which is basically the electrostatic energy of B filled by electric charge with density o in 
the presence of an ideal conductor filling the lower half-plane. As a function of the TVs it 
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obeys the infinite set of dispersionless Hirota relations. The corresponding Lax function 
performs the conformal map from the upper half plane EI to the complement of the "fat 
slit" B in H. Increasing 7\ and keeping all TVs with k > 2 fixed, one obtains a growth 
problem of Laplacian type in the upper half plane, with a specific boundary condition on 
the real line, which is associated with the dispersionless KP hierarchy in the same way 
as the problem in the whole plane is associated with the dispersionless 2DTL hierarchy. 

The tau- function (14. 5[) regarded as a function of "slow times" T k = fit k , k G Z, has 
a similar % — > limit yielding tau-function of the dispersionless 2DTL hierarchy. It 
generates a growth problem for "fat slits" in the exterior of the unit disk (Fig. [2b), which 
will be discussed elsewhere. 
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